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Abstract 

Jellium correction on the fissionability parameter is estimated within 
a Liquid Drop Model of the charged metallic cluster. This correction 
modifies the critical condition of fission and, it becomes relevant for 
small multicharged clusters. 

PACS numbers: 36.40.Qv, 36.40.-c 

Adapting the Liquid Drop Model (LDM) from nuclear physics to the descrip- 
tion of finite systems (metallic clusters) |], |2| or even infinite ones (metals) 
[||, |H is straight-forward. The physical insight of the model is explained in 
the following way || : In a system with positive surface tension and positive 
curvature energy, the energy minimizes when the system has an environment 
that is as bulklike as possible, i.e., when the surface is as small in area and 
as concave as it can be for a given volume. Thus, the total energy of an 
extended system of volume V and surface area S can be readily expressed as 
a sum of volume, surface, electrostatic (for charged systems) and curvature 
terms. 

Despite the errors arising from the shell structure of finite systems, LDM 
can be successfully used in understanding the fission of multiply charged 
metallic clusters ||. Neglecting the curvature energy (with supposition on 



a small curvature of the droplet), the total energy of the charged droplet is 
assumed to be comprised of three terms in this approach, 



E T = E b + E s + E c . (1) 

The leading term of (1) may be written = e ■ N, where iV is the number 
of constituent particles and e is the energy per particle in the limit ^ — ► oo. 
When a spherical droplet of radius R Q is deformed (conserving the volume) 
toward an ellipsoidal shape its area changes to 

S « AnR 2 (l + , (2) 

upon the second-order approximation in powers of the deformation parameter 

The surface energy, E s , in this the simplest deformation of the droplet is 
therefore 

E s « AnaR 2 (l + , (3) 

where a is the surface tension. Also, in the simplest manner, when the 
deformed droplet electrified to a charge q, it has an electrostatic (Coulomb) 
energy 

Ec « a -^— ( 1 - -P 2 ) . (4) 

Here, a = | for a uniform distribution of charge q inside the droplet or a = | 
when the charge g is spread over the surface only. Looking at eqs. (3 - 4) 
one can observe that the relative strengths of the Coulomb and surface forces 
will determine the stability of the spherical shape and hence its tendency to 
fission. 

The Coulomb energy of a charged metallic cluster is a long - standing 
controversial problem in cluster physics. Thus, the dimensionless cvoefficient 
a (see eq. 4) takes experimental values around 0.4 while from the classical 
electrostatic reasons, two different theoretical results for this coefficient have 
been derived: a = ~ from the spherical- capacitor approach |J, and a = | 
from the image-potential approach ||. Recently, Seidl and Perdew [10] have 



established that, classically, a = ~ by either approach and the experimental 
deviation from this value must be accounted for by quantum effects. Their 
theoretical result is derived in a "point-charge-plus-continuum" model for 
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jellium within which the excess electric charge is uniformly distributed in a 
surface shell between the radii (R^ — H? ) 3 and R Q . In this model, a — \ is 
a correct result in the limit z • < 1 which means large clusters and not 
so many electric charges (z ■ e). Actually, according with the above model, 
a is larger than | depending on the charge z ■ e and the number of the 
constituent particles. This follows from the fact that, increasing the cluster 
charge or decreasing the number of the cluster constituents, the volume of 
the charged surface shell increases relative to that of all the cluster volume 
and, consequently, a — >■ | which is the right value for uniformly charged 
sphere. This assumption can be relevant for small multicharged clusters. 
With these in mind somewhat corrections on the critical condition of cluster 
fission within a Liquid Drop Model (LDM) are necessary. This is, briefly, the 
aim for that the present paper is addressed. 

In the " Wigner - crystal" model of jellium, the positive charge of the ions 
is replaced by a rigid uniform positive background with the charge density 
+en = inside the sphere of radius R Q . The valence electrons are point 
charges sitting at equilibrium position inside the sphere of radius R Q . Some- 
times, by various reasons, the cluster ionization, the above model is replaced 
by the "point - charge - plus - continuum" model of jellium. In this approach, 
the electrons to be removed remain point charges while, the other electrons 
are replaced by a continuous, mobile fluid of negative charge. This system 
achieves equilibrium when the net electric field vanishes wherever the density 
of continuum negative charge does not. In this way, when the point electrons 
are still located inside the droplet, the density of continuum negative charge 
must neutralize the positive background everywhere but, inside the spheres 
of radius r s around the point electrons the density of continuum negative 
charge must vanish. As the electrons are pulled out through the surface of 
the cluster their " classical holes" deform, remaining behind inside the surface 
while conserving its volume, after the electrons have been removed to x = oo, 
these " classical holes" are spread over a shell of charge density +en between 
the radii R Q and R = (R^ — zrl)~\ where z is the number of the electrons 
which have been removed. In this way, the Coulomb potential, subject of 
the Laplace equation both inside the neutral part of the cluster (0 < r < R) 
and outside the cluster (r > R Q ) and of the Poisson equation in the cluster 
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region between R Q and R, reads 

A, if r < R; 

$(r) = { -f r 2 - f + C, \iR<r<R o] (5) 

if r>R a . 

Imposing appropriate boundary conditions, the A,B,C and D constants are 
readily determined and, the electrostatic potential in the charged shell (R < 
r < R Q ) is 

en 9 R? en / n ,\ en / , ,\ , , 
$ shell = r 2 en + — — [R 3 + 2R 3 ) — ( R 3 - R 3 ) , 6 

where e Q is the electric permittivity of free space. Under the assumption that 
e = e a , the electrostatic self-energy reads 

W = ^JdVen-$ shell (7) 

As we can see, for z = 1, the above amount is the electrostatic self-energy 
obtained by Seidl and Perdew ( equation 3 in |10| ). Further, we will note 

3 

by n and, in the second-order approximation in powers of this parameter, 
the above equation becomes 

W = ^. 1 -(l + V - + ^ + ..), (8) 

4yre.Ro 2 V 9 27 J' K ' 

which means that the value of the a parameter having corrections from the 
jellium approach of the metallic cluster reads 

a= 1 -(l + V -+ V - + ..) (9) 
2 V 9 27 J K ' 

and, it is easy to see that when n —>■ (the charged shell tends to an infinites- 
imal thin layer) then a | and, when n — > 1 (the charged surface shell has 
an increasing trend by various reasons: small clusters and/or multicharged 
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clusters) then a |. On the other hand, the radius of the spherical cluster 
can be expressed as R Q = r s N^ which means that rj = and a becomes 

(io) 

which shows that the dimensionless parameter a is larger than \ depending 
on the charge z and the number of cluster constituents N. 

As we have said, the above correction on the value of the dimensionless 
parameter a may have somewhat relevance on the critical condition of cluster 
fission within the LDM. 

For a droplet N z+ , where N is the number of atoms and z is the charge 

2 

in units e, the surface energy is proportional to iV» and is a minimum when 
the droplet is spherical. The surface energy disfavors deformation (see eq. 
3). On the other hand, the Coulomb energy depends on N~s and is a maxi- 
mum when the droplet is undeformed (see eq. 4); it favors the deformation. 
The relative strengths of the Coulomb and surface forces will determine the 
stability of the shape and hence its tendency to fission. Looking at equations 
(3 - 4) and, taking into account the last equation for the a, one can say 
that the droplet is unstable to (classical) spontaneous decay, even at zero 
temperature, when 

77( 1 + 9^ + 2^ + '"M^ ( 1 + ^ + 2^ + "'))^ n) 

647r 2 erfcr 
= e 2 ' 

Looking at the above equations we may say that the corresponding critical 
size, N c , for a given by (10) slightly differs from that corresponding to a — \. 
The jellium correction rises the lower bound of the critical sizes. 

With high enough excitation, fission is always an allowed decay channel. 
Whether fission is observed or not, however, depends on where the threshold 
energy for competing reactions lies in relation to the fission threshold. Evap- 
oration is the main competitor in both nuclear and cluster fission. Looking, 
again, at the above equations one may see that it is expected a slight de- 
creasing of the fission barrier for a given by (10) relatively to that for a=\. 

As a concluding remark, one may say that the above jellium correction on 
the fissionability criterion within a LDM, apart from the fact that promotes 
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a correct (natural) way to account for the Coulomb energy as a part in 
the fission condition, this rises (indeed not very much) the lower bound of 
the critical sizes as a consequence of a increasing and, on the other hand, 
decreases the fission barrier by the same considerations. 
The author thanks to Professor M. Apostol for many useful discussions. 
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